In number theory, the Divisor summatory function is a function that is a sum over the divisor function. It frequently occurs in the study of the asymptotic behaviour of the Riemann zeta function. The various studies of the behaviour of the divisor function are sometimes called divisor problems.
Contents |
The divisor summatory function is defined as
where
is the divisor function. The divisor function counts the number of ways that the integer n can be written as a product of two integers. More generally, one defines
where dk(n) counts the number of ways that n can be written as a product of k numbers. This quantity can be visualized as the count of the number of lattice points fenced off by a hyperbolic surface in k dimensions. Thus, for k=2, D(x)=D2(x) counts the number of points on a square lattice bounded on the left by the vertical-axis, on the bottom by the horizontal-axis, and to the upper-right by the hyperbola jk = x. Roughly, this shape may be envisioned as a hyperbolic simplex. This allows us to provide an alternative expression for D(x), and a simple way to compute it in time:
If the hyperbola in this context is replaced by a circle then determining the value of the resulting function is known as the Gauss circle problem.
Finding a closed form for this summed expression seems to be beyond the techniques available, but it is possible to give approximations. The leading behaviour of the series is not difficult to obtain. Dirichlet demonstrated that
where is the Euler-Mascheroni constant, and the non-leading term is
Here, denotes Big-O notation. The Dirichlet divisor problem, precisely stated, is to find the infimum of all values for which
holds true, for any . As of 2006[update], this problem remains unsolved. Progress has been slow. Many of the same methods work for this problem and for Gauss's circle problem, another lattice-point counting problem. Section F1 of Unsolved Problems in Number Theory [1] surveys what is known and not known about these problems.
So, the true value of lies somewhere between 1/4 and 131/416; it is widely conjectured to be exactly 1/4. Direct evaluation of lends credence to this conjecture, since appears to be approximately normally distributed with the standard deviation of 1 for x up to at least 1016.
In the generalized case, one has
where is a polynomial of degree . Using simple estimates, it is readily shown that
for integer . As in the case, the infimum of the bound is not known. Defining the order as the smallest value for which
holds, for any , one has the following results:
Both portions may be expressed as Mellin transforms:
for . Here, is the Riemann zeta function. Similarly, one has
with . The leading term of is obtained by shifting the contour past the double pole at : the leading term is just the residue, by Cauchy's integral formula. In general, one has
and likewise for , for .